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ABSTRACT 

We present new Hubble Space Telescope (HST) observations of NGC 300 taken as part of the ACS 
Nearby Galaxy Survey Treasury (ANGST). Individual stars are resolved in these images down to 
an absolute magnitude of M F814W / = 1.0 (below the red clump). We determine the star formation 
history of the galaxy in 6 radial bins by comparing our observed color-magnitude diagrams (CMDs) 
with synthetic CMDs based on theoretical isochrones. We find that the stellar disk out to 5.4 kpc 
is primarily old, in contrast with the outwardly similar galaxy M33. We determine the scale length 
as a function of age and find evidence for inside-out growth of the stellar disk: the scale length has 
increased from 1.1±0.1 kpc 10 Gyr ago to 1.3±0.1 kpc at present, indicating a buildup in the fraction 
of young stars at larger radii. As the scale length of M33 has recently been shown to have increased 
much more dramatically with time, our results demonstrate that two galaxies with similar sizes and 
morphologies can have very different histories. With an iV-body simulation of a galaxy designed to be 
similar to NGC 300, we determine that the effects of radial migration should be minimal. We trace 
the metallicity gradient as a function of time and find a present day metallicity gradient consistent 
with that seen in previous studies. Consistent results are obtained from archival images covering the 
same radial extent but differing in placement and filter combination. 

Subject headings: galaxies: evolution, galaxies: individual (NGC 300), galaxies: spiral, galaxies: stellar 
content 



1. INTRODUCTION 

The standard model of galaxy formation has the in- 
ner parts of galaxies forming before the outer parts, as 
a res ult of increas i ng timescales for gas infall with ra- 
dius (lLarsonlll976t iWhite fe Frenkl 119911 : iBurkert et al.1 
[1991 IMo et al.1 119981: INaab fe Ostrikerl 12006^ This 
"inside-out" growth scenario has likewise been seen 
in V-body/Smooth Particl e Hydrodynamics s imulations 
of disk galaxy evolution (|Brook et al.l H006), and has 
been used to explain the observed abundance gradi- 
ents ij2__the_^lfikY__Wayvi^ L _che mod- 
els dMatteucci fe Francoii 119891: iChiappini et al.l Il997t 
iBoissier fe Prantzoslll999T ). 

There are multiple mechanisms which could contribute 
to an observed trend of star formation taking place more 
recently in the outer parts of galaxies. The first, as pre- 
viously mentioned, is that gas does not accumulate in 
the outer disk until later times, or that it accumulates 
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more slowly such that the total mass density of the outer 
disk increases with time. Alternatively, if the gas is in 
place but the star formation timescale in the outer disk 
is longer, the stellar mass density would increase more 
slowly with time in the outer disk. 

If disks do form "inside-out," one would expect to see 
negative radial gradients in age and metallicity. Age gra- 
dients would result from a higher percentage of older 
stars in the centers of galaxies, where star formation 
started earlier. This older population of stars would en- 
rich the surrounding gas, causing more recently formed 
stars toward the center to be more metal-rich than stars 
formed in the outskirts of the disk. 

Negative gradients have been seen in several surveys 
of nearby galaxies, and they appear to be common, if 
not universal, among larg er disks. In the sample o f 
nearby galaxies presented in lMuhoz-Mateos et al.l (2007) , 
FUV — K color was used as a proxy for star formation 
rate (SFR), since it represents the ratio of young to old 
stars and shows what fraction of the total star forma- 
tion has been recent. Large spirals were found to have 
predominantly negative gradients, while low-luminosity 
systems show a considerable scatter in the slopes of their 
gradients, with both positive and negative values, in- 
dica ting a wider range of p ossible formation histories. 
The iMacArthur et" al. (2009) study of nearby spirals us- 
ing spectral synthesis finds both negative gradients and 
flat profiles in stellar age and metallicity. 

A flat radial profile in age and metallicity may be con- 
sistent with inside-out growth if gradients established 
early were erased by interactions or subsequent radial 
migration. Especially in high-mass systems, stars can 
scatter off of spiral struct ure and change their radii while 
retaining circular orbits (|Sellwood fc Binnevl l2002'). The 
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TABLE 1 

Summary of NGC 300 and M33 properties 





NGC 300 


M33 


Distance 


2.0 Mpc a 


800 kpc b 


Type 


SA(s)d c 


SA(s)cd c 


M B 


-17.66 a 


-18.4 d 


Scale length (K) 


1.3 c 


1.4 e 


Circular velocity 


97 km s _1 f 


130 km s" 1 s 


Estimated mass 


2.4 X 10 10 M Q f 


5 x 10 10 M s 



a IDalcanton et alj j2009T l 
b IWilliams et al.l Il2009al) 

c NED 

d IVila-Costas fc Edmundl lfl992 
e IMunoz-Mateos et al.l (120071 ) 
f Puchc et al. (1990) 
8 Corbclli & Salucci (2000) 



simulations of Roskar ct al. (2008) showed that some per- 
centage of older stars currently residing in the outer disk 
of a galaxy actually formed closer to the center of the disk 
and migrated outward to their current location. Thus, 
an observed older stellar population in the outer disk of 
a massive galaxy may not necessarily indicate that the 
outer disk formed early. 

The ideal galaxies for studying the evolution of disks in 
the absence of major mergers are undisturbed, pure-disk 
(i.e., bulgeless) systems that maintain weak spiral struc- 
ture to suppress radial migration. While one method 
of studying galaxy evolution is to compare high-redshift 
galaxies with their local counterparts, difficulties with 
this method include the dramatic falloff of surface bright- 
ness with redshift and the inability to trace the history 
of a single galaxy. Color-magnitude diagrams (CMDs) 
derived from resolved stellar populations can provide de- 
tailed, time-resolved studies of the disk evolution. This 
method restricts the possible targets to relatively local 
galaxies. 

NGC 300, in the Sculptor Group, is the nearest iso- 
lated late-type disk galaxy, and is thus an ideal target. 
IMunoz-Mateos et~a l. (2007) have already suggested an 
inside-out growth scenario for NGC 300 based on its 
broadband colors. The resolution of the Hubble Space 
Telescope (HST) enables us to use CMD fitting to find 
the star formation histories (SFHs) of galaxies outside 
the Local Grou p. The ACS Nearby Ga laxy Survey Trea- 
sury (ANGST; IDalcanton et al.l I2009D was designed to 
create a volume- limited sample of nearby galaxies, allow- 
ing for an unbiased accounting of star formation in the 
nearby universe. The continuous radial strip of NGC 300 
imaged as part of ANGST allows us to recover the SFH 
of this galaxy out to ~ 5 kpc. 

Evidence for inside-out growth (a decrease in mean 
stellar age with radius) has also been found in M33, 
a galaxy very similar to NGC 300 (see Table Q] for 
a comparison of thei r main properties). In M33, 
IWilliams et al.l (|2009aD and Holtzman et al. (in prep.) 
used resolved stellar populations from four HST /ACS 
fields to derive SFHs using CMD fitting. From these 
SFHs, they infer the stellar surface density of the disk at 
different times throughout the galaxy's history and find 
significant evolution in the scale length of the disk. The 
increase in scale length with time, indicating that the 
SFR has been increasing in the outer disk, is suggestive 
of inside-out growth. 




Fig. 1.— DSS image of NGC 300 with ACS fields overlaid. 
The three contiguous ACS fields forming a radial strip are from 
ANGST, while the other fields are archival data. 

While outwardly similar, NGC 300 and M33 may 
have different histo ries. M33 ha s a d isk break at 
~ 6 scale lengths (jFerguson et al.1 120071). which is a 
comm on feature in spiral galaxies dPohlen fc Truiillol 
I2006T) . However , iBland-Hawthorn et al.1 (|2005l ) showed 
that NGC 300 has a pure exponential disk out to ~ 10 
scale lengths. There are environmental differences as 
well: an Hi bridge between M33 and M31 is suggestive 
of a history of interact i on between these two galaxies 
(jBraun fc Thilkerl l2004t iBekkil |2008 f) . a finding that is 
confirmed by the distribution of red giant branch (RGB) 
stars surrounding the two galaxies (IMcConnachi e et al.l 
I2009D and the evidence fo r tidal disruption of M33's gas 
disk (jPutman et al.ll2009D . In contrast, NGC 300 is fairly 
isolated on the Sculptor filament, wit h only dwarf galax- 
ies nearby ()Karachentsev et al.l l2003). These differences 
may imply significantly different SFHs. 

In fj2l we describe our data and reduction; in $3l we 
describe our methods for determining the star forma- 
tion history and present our results; we discuss their 
implications for disk growth along with possible caveats 
and compare our metallicity results to other observations 
in 21 and we conclude with JJSj Archival data is pre- 
sented in Appendix |A"1 We adopt a WMAP cosmology 
(jSpergel et al.|[2007l) for all conversions between time and 
redshift. 

2. DATA AND PHOTOMETRY 
2.1. ACS Imaging 

As part of ANGST (GO-10915 ) we took HST Ad- 
vanced Camera for Surveys (ACS: IFord eTall fl998l) ob- 
servations of NGC 300 during 2006 November 8-10. 
Three slightly overlapping fields were observed along a 
radial strip from the center of the galaxy into the disk. 
Each field was observed for 1488s in FA75W, 1515s in 
F606W, and 1542s in .F814W. An additional deep outer 
field was planned as part of the ANGST program, which 
would have given more extended radial coverage of the 
disk; unfortunately, ACS failed before this observation 
was obtained. 

An additional six archival ACS images, scattered 
across the disk of NGC 300, are presented in Appendix lAl 
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Fig. 2.— Color image of the radial strip of NGC 300 observed as part of ANGST (red: F814W, green: F606VK, blue: F475W). 
Inclination-corrected radial bins of width 0.9 kpc are overlaid as white lines, with radii in kpc labeled. 



We include the archival data as a useful consistency check 
that the ANGST data is representative of the disk as a 
whole. However, since the observing conditions were not 
identical (different exposure times, different filters, and 
no spatial continuity) , we keep the analysis of the archival 
data separate from that of the new observations. These 
images were originally taken as part of the Araucaria 
project to determine Cepheid distances to nearby galax- 
ies (GO-9492), so the fields were selected to sample the 
Cepheid population at different galactocentric distances 
and ar e therefore placed on more active star- forming re- 
gions (|Bresolin et all [200l[RTzzi et al.|[200l) . The ob- 
servations were obtained between the dates of 2002 July 
17 and December 25, with exposure times of 1080s in 
FAZbW and F555W, and 1440s in F8UW. Each ob- 
servation was split between two exposures for cosmic-ray 
removal and coverage of the chip gap. The exposures 
were calibrated and flat-fielded using the standard HST 
pipeline. The locations of all ACS fields are shown in 
Figure [TJ 

2.2. Photometry 

While exposures were combined using MultiDrizzle to 
produce images in this paper, photometry was done us- 
ing all individual exposures simultaneously. For photom- 
etry, we use DO LPHOT, a modified version of HSTphot 
(Dolphin 2000) optimized for the ACS. DOLPHOT fits 
the ACS point spread function (PSF) to all of the stars in 
each exposure, determines the aperture correction from 
the most isolated stars, combines the results from all ex- 
posures, and converts the count rates to the Vega magni- 
tude system. Details of the photometry and quality cuts 
us ed for the ANG S T samp le an d archival data a r e give n 
in IWilliams etHI (|2009bl i and IDalcanton et all (2009). 
We require that stars in the final sample are classified as 
stars, not flagged as unusable (too many bad or saturated 
pixels or extending too far off the edge of the chip) , have 
S/N > 4, and have \sharppQQe\y + sharppsuwl < 0.274. 
Sharpness indicates whether a star is too sharp (per- 
haps a cosmic ray) or too broad, and these cuts exclude 
non-stellar objects (such as background galaxies) that 
escaped the earlier cuts. We also cut on the crowding 
parameter, which is defined as how much brighter a star 
would have been measured if nearby stars had not been 
fit simultaneously. Stars with a high crowding parame- 



ter are more likely to have erroneous photometry, but a 
very strict cut has the effect of preferentially removing 
young stars, since these are usually found in clusters. We 
require crowdFw&w + crowdFsuw < 0.6 mag. 

We characterize the completeness of our sample in 
terms of magnitude, color, and position by inserting at 
least ~ 10 6 artificial stars (~ 50, 000 stars at a time) in 
each ACS field. DOLPHOT is also used to perform the 
artificial star tests. Individual stars are inserted into the 
original images and their photometry is rc-measured. Ar- 
tificial stars are labeled as "detected" if they were found 
by DOLPHOT and met the quality cuts described above. 
The 50% completeness limit for ANGST ranges from 
F606W = 26.7, F8UW = 25.8 in the crowded center of 
the galaxy to F606W 7 = 28.4, F8UW = 27 A in the outer 
regions. For the archival data, 50% completeness limits 
range from F555W = 26.8 to 27.7 and F8UW = 25.8 
to 27.2. 

2.3. Identification of Duplicate Stars 

The ANGST fields overlap by - 7" so that the fields 
could be aligned into a single radial strip. When com- 
bining the photometry of all stars into a single catalog, 
we identify and remove duplicate stars in overlapping re- 
gions by binning star positions in each field and find- 
ing bins that contain stars from more than one field. 
The bin size is set at 1", big enough so that nearly all 
bins in the overlap regions contain stars in both fields 
but small enough so that the edges of the overlap re- 
gions are fairly smooth. The three ANGST fields are 
designated WIDE1, WIDE2, and WIDE3, with WIDE3 
containing the center of the galaxy and WIDE1 farthest 
from the center. For the WIDE1-WIDE2 overlap, stars 
in WIDE1 are kept, and for the WIDE2-WIDE3 overlap, 
stars in WIDE2 are kept. The World Coordinate System 
(WCS) for each of 3 fields is provided by the standard 
HST pipeline, but there are very small offsets between 
fields. We found that relative to WIDE2, WIDE1 has an 
offset of Aa = 0.23", AS = 0.036", and WIDE3 has an 
offset of Aa = -0.13", AS = -0.14". 

2.4. Dividing Stars into Radial Bins 

To assign an inclination-corrected galactocentric dis- 
tance to each star, we assume the following galaxy pa- 
rameters: a = 13.722833°, 5 = -37.684389° (galaxy 
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Fig. 3. — CMDs of radial bins from ANGST data. Only stars with S/N > 4 and high-quality photometry have been included (see JJ22J. 



cente r), i = 42° (inc lination), 6 = 111° (position an- 
gle) JK mi etHll200l . To convert radius on the sky to 
physical distance, we assume a distance of 2.0 Mpc to 
NGC 300, which is the distance we measure using the ti p 
of the red giant branch (TRGB) (jDalcanton et al.ll2009j) . 
The maximum radius we find for a star in the ANGST 
fields is r = 5.35 kpc. 

Stars are divided into 6 radial bins based on their 
(inclination-corrected) distance from the center of the 
galaxy. Artificial stars are divided in the same way. Our 
bins are spaced 0.9 kpc apart in radius, to encompass all 
of the observed stars in annuli of equal width (Figure [2]). 
The spacing of the bins is motivated by the length over 
which we expect ra dial mixing to blur the populations 
(Rosk ar" et al.l l2008). The observed area is calculated for 
each bin so that we can derive the correct surface density 
of star formation. CMDs for stars in each bin are shown 
for the ANGST data in Figure El 

In the various fits to distance shown in $4] and Ap- 
pendix [A] the distance used for each radial bin is the 
mean distance from the center of pixels in that region. 

3. STAR FORMATION HISTORY ANALYSIS 
3.1. Method 

To derive the SFH of each radial bin, we use the 
well-established technique of comp aring the observe d 
CMP to a set of mode l CMDs (e.g.. iGallart et allll99" 
Hernandez et al.l 119991: iHoltzman et all Il999t fpolphi 



2001 ISkillman et all [20031: Harris k Zaritekvl [200l 
Gallart et al.ll2005[ ). Typical fitting codes use stellar evo- 



lution models that predict the properties of stars of dif- 
ferent masses for a range of ages and metallicities. From 
the predicted luminosity and temperature, the magni- 
tudes of the stars are determined for a given filter set. 
Stars are then placed on a synthetic CMD following the 
mass distribution of an assumed initial mass function 
(IMF) and binary fraction for each age and metallicity. 
With distance and extinction either set or included as 
additional free parameters, these model CMDs are lin- 
early combined until the best fit to the observed CMD is 
found. The ages and metallicities of the CMDs that went 
into the best fit tell us the ages and metallicities of the 
underlying stellar population, while the weights given to 
the CMDs provide the SFR at e ach age. 

We use MATCH, described in iDolphinl (pOOl . to de- 
rive the SFH for each radial bin. This code finds the 
maximum-likelihood fit to the CMD assuming Poisson- 
sampled data. We assume an IMF with a slope of -2.35 
(|Salpeterl 119551 ) between 0.1 and 120 M© and a binary 
fraction of 0.35. Given that our CMD only includes stars 
with masses > 1 M©, adopting a single Salpeter slope is 
likely to be a valid assumption. While a power-law IMF 
was used by MATCH for comp utational ease , the SFRs 
we show have been scaled to a lKroupal (|200lD IMF. This 
is possible because, for stars massive enough to be ob- 
served, the two IMFs are very similar. The choice of 
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Fig. 4. — Lifetime SFHs of radial bins from ANGST data. Error bars are the quadrature sum of th e uncer tainties from distance and 
extinction and the 68% confidence interval from Monte Carlo simulations. SFRs have been scaled to a[Kroupa (2002) IMF. The scale of 
the vertical axis differs in each panel. 



IMF affects only the normalization of the SFH, but not 
the time dependence, so relative differences among age 
and radial bins do not depend on the IMF. 

Synthetic CM Ds are constr u cted from the theoreti- 
cal iso chrones of iGirardi et al.l (|2002t ) and iMarigo et al.l 
(|2008f) for ages in the range 4 Myr-14 Gyr. The 
isoch rones younger than ~ 6 x 10 7 yr were adopted 
from IBertelli et al.l (119941, . with t ransfo rmations to the 
ACS system from IGirardi et al.l (|2008l ). Age bins are 
spaced logarithmically, since the CMD changes much 
more rapidly at young ages than at old ages. Metallicity 
is allowed to vary in the range —2.3 < [M/H] < 0.1, but 
is not allowed to decrease with time. See N4.3I for a fur- 
ther discussion of metallicity and uncertainties therein. 

MATCH determines the best fit for distance and ex- 
tinction by testing a range of values. The range for 
the distance modulus is 26.3 < m — M < 26.7, span- 
ning the range of values report ed in the literature (e.g. , 
Butler et aj] 120041 ISakai et al.l[200l I G ieren et al.l 120051: 



Rizzi et al.l 120061 : iDalcanton et all 120091 ) and the extinc- 
tion range is 0.05 < Ay < 0.5. Additionally, up to 0.5 
mag of differential extinction is applied to young stars 
(< 100 Myr), since these stars are mo re likely to be 
found in dusty st ar form ation regions (iZaritskvl Il99l 
iZaritskv et alil2002[ ). The lSchlegel etall (pgl Tralue for 
Galactic extinction is Ay = 0.042 in the line of sight to 
NGC 300, but we expect the total value to be higher due 
to local extinction within NGC 300 itself. Additionally, 
we expect that dust content may vary across the extent 
of the disk, so extinction may be different in each radial 
bin. 

Completeness and observational errors are accounted 
for by including the results of artificial star tests. We 
identify the artificial stars that were placed within each 
radial bin and supply MATCH with their input and out- 



put magnitudes and whether they were detected above 
the quality cuts of our photometry. The density distri- 
bution of artificial stars mirrors the density of detected 
stars within each bin, so MATCH accounts for any radial 
variation in crowding. 

We use the V + I equivalent filter sets (F606W + 
F8UW, F555W + F8UW) to derive the SFHs in this 
paper, due to the greater depth of the F606H^ and 
F555W data as compared to the F457W and F 435VF 
data. As a consistency check, we also derived the SFHs 
using F475PF + F814W/ for the ANGST fields and found 
SFHs that were consistent within the error bars (which 
are larger in FA75W + F814W, especially in the inner 
regions) , so our color choice does not appear to affect our 
conclusions significantly. 

We assess uncertainties due to Poisson sampling of 
underpopulated regions in the CMD by running Monte 
Carlo simulations as follows: for each region, we sample 
stars at random from the observed CMDs until we reach 
the same number of stars as observed. These stars are 
then given as the input to MATCH with the distance 
and extinction fixed, and the resulting SFH is compared 
to the SFH from the original data. We repeat this pro- 
cess 100 times and define our sampling error as the val- 
ues which encompass 68% (lcr) of the Monte Carlo tests. 
Our final error bars are the quadrature sum of this er- 
ror and the systematic errors from fitting the distance 
and extinction. Our error bars do not include systematic 
uncertainties in the stellar evolution models. 

3.2. Distance, Extinction, and Crowding Effects 

Because our data span a large range of stellar densi- 
ties, the effective completeness varies strongly with ra- 
dius. To deal with this variation, we pursued two inde- 
pendent methods. The first approach is to consider only 
the portion of the CMD that is 50% complete for each 
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Fig. 5 . — Re cent (< 1 Gyr) SFHs of radial bins from ANGST data. Errors are calculated as in Figure 3] SFRs have been scaled to a 
Kroupa (2002) IMF. The scale of the vertical axis differs in each panel. The bin at 2.7 < r < 3.6 kpc contains a prominent spiral arm. 



radial bin, so that larger portions of the CMD can be 
used in the outer, less crowded regions. The second ap- 
proach is to find the completeness in the most crowded 
region (the center of the galaxy) and apply that magni- 
tude cut to every region uniformly. This removes useful 
data in the outer regions, but ensures that variations we 
see in the SFH across the disk are not due to variations 
in photometric depth. We have employed both methods 
and compared the results. 

When we allow the completeness cut to vary with ra- 
dius, the photometry in the outer four radial bins is com- 
plete to well below the red clump (see Figure [3]). This 
depth allows for a more consistent distance measurement 
across these four fields; derived distance modulus values 
for the ANGST data in the six radial bins were as follows: 
m-M = 26.61±0.07,26.60±0.07,26.43±0.09,26.43± 
0.09, 26.42 ±0.09, 26.43 ±0.09. All of the four outer bins, 
with depths that resolve the red clump, produce a con- 
sistent estimate for the distance of m — M = 26.43 ±0.09. 
The inner fields, where crowding limits the depth, have 
a discrepant distance. We hereafter fix the distance to 
m — M = 26.43 in all fields. Even if this distance is not 
absolutely correct, what is important is that setting the 
distance to the same value for all fields ensures that any 
relative changes in the positions of stars on the CMDs 
between radial bins are interpreted as differences in SFR 
and/or extinction. 

Comparing the SFH when depth is allowed to increase 
with radius and the SFH where depth is at a fixed mag- 
nitude throughout, the recent SFHs for ages < 1 Gyr 
are very similar, probably because the recent SFH is 
constrained by younger, brighter stars and is not sig- 
nificantly affected by the elimination of stars near the 
bottom of the CMD. For ages > 1 Gyr, the most sig- 
nificant difference is a 40-80% increase in the SFR at 
8-14 Gyr for the outer 3 radial bins when restricted to a 



uniform shallow depth. There is a corresponding overall 
decrease in the more recent SFR (1-8 Gyr). These tests 
indicate that using these time bins, the central region's 
SFH may be more biased to older ages than would be 
measured with less crowded data. 

To address this effect, we determined bin widths 
that reflect o u r sens itivity to age as described in 
iWilliams et~aT1 ((2009c). The resulting time bins are 
wider at intermediate ages (1-10 Gyr) for the most 
crowded regions, and thus robust against uncertainties 
in SFR within this age range. Since the most recent 
time bins contain small numbers of stars, assessing the 
statistical significance of these bins is more difficult, and 
thus we present the past ~ 80 Myr as a single time bin 
in all regions. Changes in SFR on shorter timescales do 
not affect any of the conclusions of this paper. 

Given that using all available data allows us to better 
constrain the SFH in the outer regions, the SFHs pre- 
sented in this paper are those derived allowing the 50% 
completeness limit to vary with radius. 

The mean extinction values found for the ANGST data 
are Ay = 0.10 ± 0.05 for all bins. Note that these val- 
ues do not include the 0.5 mag of differential extinction 
applied only to young stars (< 100 Myr). We discuss 
additional tests applying differential extinction to older 
stars as well in H4.2.2I 

3.3. Radially Resolved Star Formation Histories 

The SFH for each radial bin, as derived by MATCH, 
is shown in Figures 2] and O the former shows the life- 
time SFH, and the latter focuses on the recent SFH (< 1 
Gyr). The overall behavior of SFR vs. age changes with 
radius; central regions have a higher percentage of old 
stars, and the percentage of young stars increases with 
radius. The bulk of a spiral arm, visible on the color im- 
age in Figure [2l is within the radial bin at 2.7 < r < 3.6 
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Fig. 6. — Cumulative star formation (i.e., fraction of stars formed 
vs. age of galaxy) for all radial bins for ANGST data. Bin edges 
have been offset to avoid overlapping error bars. 

kpc. This region also has the most dramatic increase of 
SFR in the most recent time bin as compared with its 
level in the past Gyr, as can be seen in Figure [SI 

Appendix [S] presents the SFH for archival data in iden- 
tical radial bins. Rather than a continuous radial strip, 
these fields are scattered throughout the disk. The gen- 
eral agreement between the radial SFHs of fields selected 
in two different ways illustrates that we are not biasing 
the results by looking at only a portion of the entire disk. 

We present the full SFH for all radial bins, including 
metallicities, in Table [2] Metallicity will be discussed 
further in gOj 

4. DISCUSSION 
4.1. Stellar Disk Growth 

We can use the SFH to infer the past structure of the 
disk. There is evidence for "inside-out" growth (Fig- 
ure S]), as early star formation was more prominent in 
the inner disk. These trends are quantified in Figure [6l 
where we plot the cumulative SFH. Cumulative plots are 
presented at the full resolution of the CMD fit, since 
changes in SFR due to uncertainties typically occur in 
adjacent time bins and thus have a very small effect on 
the cumulative SFH. The outer parts of the disk formed 
a greater fraction of their stars at recent times than the 
inner parts of the disk, consistent with the "inside-out" 
growth scenario. However, even the outer regions are 
fairly old, with ~ 75% of stars formed by 4 Gyr ago. 

We emphasize that when we refer to our "outer" disk 
observations in this paper, the galactocentric distances 
are still within 5.4 kpc for a gala xy that has been shown 
to ex tend to at least 14 kpc IjBland-H awthor n et al.l 
2005). Therefore, the reader should keep in mind that 
all distinctions between "inner" and "outer" we can make 
with our data are still within what some might consider 
the "inner" regions of NGC 300. 

We can also infer past stellar surface density if we as- 
sume that the majority of stars formed within the radial 
bin in which they are currently found (see t j4.2.3[) . We 
calculate the stellar mass in each radial bin that formed 
before a given time by summing all the mass formed up 
to that time in the derived SFH. For the purposes of 
calculating stellar surface density we require a uniform 
set of time bins for all observed regions, so we choose a 
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Fig. 7. — Stellar surface density versus radius for multiple time 
bins. Points have been offset from one another by 0.05 kpc to avoid 
overlapping error bars. T o compare this surface density with one 
assuming a Kroupa (2001) IMF, one should divide by 1.5. 

coarse binning scheme of 4-100 Myr, 100 Myr-1 Gyr, 1- 
5 Gyr, 5-10 Gyr, and 10-14 Gyr. The resulting surface 
mass density profile is shown in Figure [7J We omit the 
most recent time bin from the plot, as the surface density 
has not changed significantly in the past 100 Myr. The 
stellar surface density has changed more substantially in 
the outer regions of the disk than in the central regions. 

For each time bin, we fit an exponential disk model to 
derive the scale length of the disk; the results of these 
fits are shown in Table The scale length increased 
slightly over the galaxy's history, from 1.1 ± 0.1 kpc at 
early times to 1.3 ± 0.1 kpc by 1 Gyr ago. However, 
our error bars are also consistent with no scale length 
evolution. A predominantly old disk is in contrast to 
the dramatic changes in scale length than has been seen 
in M33 or has been inferred fro m in situ studies of disk 
evolution at high redshift (e.g., rTruiillo fc Pohlenl [20051 : 
iBarden et aIll20Q5l: lAzzollini et al.ll2008H . M33 has been 
shown to have a scale length th at increases by nearly a 



factor of 2 inside the disk break (Williams et al 



in agreement with the predictions of IMo et al 



2009a), 
119981 

Although M33 is a near twin of NGC 300 in mass and 
morphology, NGC 300 is differen t in that it lacks a disk 
break in its exponential profile (Bland-Hawthorn et al. 
I2005f ) and is much more isolated than M33. 

The value for the scale length derived by sum- 
ming over the stellar mass formed in the derived 
SFH agrees remarkably well with the scale length for 
the fC-band stellar mass surface density of 1.29 j^'j^ 
(M uiioz-Mateos et al. 2007). Since massive stars have a 
high mass-to-light ratio in the -ftT-band, agreement with 
the scale length as traced by the stellar mass formed 
in the disk is expected for a galaxy dominated by old 
stars (> 5 Gyr). Scale lengths measured at shorter wave- 
lengths are predictably larger, sinc e these trace you nger 
stellar populatio ns: 1.47 kpc in I (|Kim et al.ll2004D and 
2.17 kpc in Bj (|Carignanl TT985. scaled to a distance of 
2.0 M pc). M33 has a similar scale length in the -ftT-band 
of 1.4 ijMunoz-Mateos et al.ll2007D . In the outer fields of 
M33 (4-6 kpc, or ~ 3-4 scale lengths), only ~ 20% of 
the stars formed by 8 Gyr (z ~ 1), whereas in NGC 300, 
50-70% of stars in the equivalent outer bins had formed 
by this time. Thus, the disk of NGC 300 appears to 
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Fig. 8. — Cumulative star formation for the entire galaxy out to 
5.4 kpc. Shaded region indicates uncertainty. 

be older overall than M33. If NGC 300 experienced sig- 
nificant inside-out growth, it may have happened earlier 
than we are sensitive to with CMD fitting. 

Figure E shows the cumulative SFH for all of NGC 300 
out to 5.4 kpc. For each radial bin, we assumed the ob- 
served SFH was characteristic of the entire annulus, and 
summed stellar mass formed over all radial bins. Overall, 
~ 75% of stars in this portion of NGC 300 formed by 8 
Gyr ago. 

4.2. Complicating Effects 
4.2.1. Photometric Depth 

Since the red clump is not resolved in the crowded in- 
ner regions, the central SFH (older than ~ 1 Gyr) is 
less certain (see §3.2j) . and could possibly affect our scale 
length derivation. As an alternative, we fit an exponen- 
tial profile to only the outer four radial bins, and found 
essentially the same results as the fit to all the regions, 
albeit with larger error bars. Thus the crowded inner re- 
gions are not significantly affecting the scale lengths we 
derive. 

To further investigate whether the SFH found for 
M33 is consistent with the observed stellar populations 
of NGC 30 0, we constructed mo del CMDs from the 
M33 SFHs (IWilliams et all l2009al) using the photome- 
try statistics of our NGC 300 data. For each radial bin 
in NGC 300, we used the SFH from the observed field 
in M33 that most closely matched in radius. We ran 
MATCH on the NGC 300 data while enforcing the M33 
SFH model, scaling the SFH overall so that the number 
of stars produced by the model was equivalent to the 
number of stars observed in NGC 300. In all cases, there 
were substantial residuals and an obvious mismatch be- 
tween the NGC 300 data and M33 SFH model. 

Figure |H] compares the original MATCH fit to the 
NGC 300 1.8 < r < 2.7 kpc radial bin with the results 
of fitting the M33 model for 2.5 kpc to the same data. 
The top series of panels shows the MATCH fit to the 
NGC 300 data, along with the weighted residuals (data 
minus model divided by Poisson noise). Residuals for 
the best fitting model are very small. The bottom se- 
ries of panels shows the M33 model CMD at this radius, 
along with the (very substantial) weighted residuals. We 
present this radial bin as an example; residuals were com- 



parably large in other bins. Figure 8 in William s et al.l 
(|2009bl) demonstrates the effect of age and metallicity 
on the red clump and asymptotic giant branch (AGB) 
bump: higher metallicity gives redder colors and fainter 
magnitudes, while older ages result in an even stronger 
push towards fainter magnitudes. As shown in Figure [51 
the RGB is redder in NGC 300 than in M33, and the red 
clump is fainter, consistent with an older, more metal- 
rich stellar population. This comparison confirms that 
the NGC 300 data is not consistent with an M33-likc 
SFH, but rather that NGC 300 formed the majority of 
its stars much earlier than M33. 

4.2.2. Extinction 

Given what appear to be dusty regions visible toward 
the center of NGC 300 (Figure H}, it is possible that dif- 
ferential extinction may be a factor for all ages of stars. 
Our SFH derivation accounts for differential extinction 
of young stars (< 100 Myr), but if there are dust lanes 
present, differential e xtinction could affect older stellar 
populations as well. iRoussel et al.l (|2005l ) studied ex- 
tinction in Hn regions in NGC 300 and found values 
for A(Ha) ranging from 0.15 to 1.06. While we expect 
the Hn regions to contain predominantly young stars, 
we can use their estimates as a rough guide for assessing 
extinction. In MATCH, the differential extinctions for 
young stars and for the stellar population in general are 
set by two different parameters, with the total extinc- 
tion for young stars randomly assigned up to the sum of 
these two values. Given that we set the maximum extinc- 
tion for young stars to be 0.5, if we additionally supply a 
value for all stars of 0.5, the young stars wo uld receive ex- 
tinctio n roughly in the range observed by IRoussel et al.l 

Pol . 

We experimented with adding extinction parameters 
from 0.1 to 0.5 mag, meaning that each star is randomly 
assigned an extinction from zero to an upper limit of 
the extinction parameter. We found that the fit qual- 
ity decreased with the amount of differential extinction 
added. In the central bin, where the presence of dust 
is clearly visible, this difference was minimal: up to 4%. 
In the remaining bins, the fit was up to 25-40% worse. 
For the central bin, the next best fit after no additional 
differential extinction is when this parameter is set at 
0.2 magnitudes, so we recalculated the surface density at 
each time bin using up to 0.2 mag of additional extinc- 
tion for only the central bin. Within the error bars, the 
scale length evolution remained the same. 

4.2.3. Radial Stellar Migration 

iWielenl (|1977l ) suggested that substantial migration 
among the stars i n a d isk galaxy is expected, and 
iSellwood fc Biniievl (|2002T ) showed that transient spirals 
can provide an effi cient mechanism to redistribute stars 
radially. Recently, IRoskar et al.l (|2008|) have reproduced 
this phenomenon in simulations of growing disks. There- 
fore, the radius at which stars are obser ved may not 
be th e radius at whi ch they formed (e.g ., IWielen et al.l 
I1996T) . According to IRoskar et al.l (|2008D . migration ef- 
fects become more important with increasing radius, in 
terms of the fraction of stellar mass that is composed of 
migrated vs. in-situ stars. In their Milky Way-analogue 
model, migration significantly influenced the surface den- 
sity beyond ~ 2 scale lengths. This suggests that the 
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Fig. 9.— Top panels: results of fitting model CMDs to the NGC 300 data. Left to right: NGC 300 binned CMD, best fit CMD to the 
NGC 300 data as found by MATCH, weighted residuals ((data - model)/Poisson noise). Bottom panels: results of fitting a synthetic CMD 
created from the SFH of the M33 field at ~ 2.5 kpc to the NGC 300 CMD at 1.8-2.7 kpc. Left to right: NGC 300 binned CMD, best fit 
CMD to the M33 data, weighted residuals. The grayscale is defined as (NGC 300 data) - (M33 model)/ (Poisson noise), with darker areas 
being higher values (i.e. black areas show where there are more stars in NGC 300 and white areas show where there are more stars in 
M33). NGC 300 has a redder RGB and a fainter red clump, which is consistent with an older, more metal-rich stellar population. 



three outer radial bins in NGC 300 may contain stars 
that have migrated from nearer the center of the galaxy. 
Migration of stars from the inner to the outer disk could 
potentially erase the signature of inside-out growth in 
NGC 300, because radial bins could be contaminated by 
stars that did not actually form within them. If NGC 300 
has had more substantial migration than M33, that could 
explain the discrepancy between the scale length evolu- 
tion of the two galaxies. On the other hand, both M33 
and NGC 3 00 are significantly l ess massive than the disks 
used in the lRoskar et al.l ()2008l ) simulations, so the effect 
of radial mixing may be less severe, due to weaker spiral 
structure. 

To examine the effects of migration in NGC 300, we 
ran an A^-body simulation of a disk galaxy designed to 
be similar in mass and angular momentum to NGC 300. 
The mass at 11.8 kpc (assuming a distance of 2.0 Mpc 
to NGC 300 ) was measured to be 2.4 x 10 10 M Q by 
iPuche et al.l (1990) using velocity measurements from 
Ha. However, the disk of NGC 3 00 actually extends sig- 
nifica ntly farther than this fe.g.. iBlan d-Haw thorn et al.l 
120051) . We estimate the total mass of NGC 300, includ- 
ing baryonic and non-baryonic mass, to be ^lO 11 Mq for 
the purposes of the simulation. 

To estim ate the spin parameter A, we use th e formula 
derived bv lHernandez fc Cervantes-Sodil (|2006D : 



as Rd = 1-47 kpc, an d the velocity of Vd = 93.4 ±8.2 km 
s _1 was measured bv IPuche et al.l ([1990) with Hi obser- 
vations. Putting these values in equation [TJ we find that 
for NGC 300, A « 0.0355. The Milky Way model (MW) 
uses a total mass of 10 12 Mq and a spin of A sa 0.039. 
Both models are evolved to 10 Gyr; the scale lengths of 
the two models at the end of 10 Cyr are 1.9 kpc and 3.5 
kpc for N300 and MW respectively. 

Comparing the results of the s imulated NGC 300 with 
the Milky Way-sized galaxy from Roskar ct al. (200H), we 
find that the effect of radial migration is substantially 
reduced in the smaller galaxy. In Figure I10[ we show 
the probabilities of migration for the Milky Way model 
(top) and the NGC 300 model (bottom). The left-hand 
panels should be interpreted as, "if a star is currently at 
Rfinai, what is the probability that it formed at i?/ ™?" 
The right-hand panels, conversely, show "if a star forms 
at Rform, what is the probability that it will end up 



at R 



final 



?" 



A = 21.6 



■Rd/kpc 



{V d /km s" 1 ) 3 / 2 ' 



(1) 



where Rd is the disk scale length and Vd is the rota- 
tion velocity of a flat rotation curve. The J-band scale 
length of NGC 300 was measured by Kim et al.l (|200l 



In all panels, probability is indicated by 
color, with redder colors meaning higher probabilities. 
In the MW model, most stars form at < 12 kpc (top 
right panel), but these stars are found at radii up to 16 
kpc (top left panel). Stars at all radii have a tendency 
to end up farther out than where they formed. In the 
N300 model, however, most stars are found near their 
birth radii. 

Figure [10] shows that particles in the N300 model 
undergo much less radial migration than those in the 
more massive MW model. Since migration occurs when 
stars scatter off disk asym metries such as spiral arms 
(jSellwood fe; Binnevl I2002D , weaker spiral structure re- 
duces the probability of migration. A disk must remain 
kinematically cool to sustain recurrent transient spirals 
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Fig. 10. — Probability of migration, for MW model (top) and N300 model (bottom). Left-hand panels show, for a star at a given Rfi na i, 
the probability that it formed at Rf orm . Right-hand panels show, for a star at a given Rf orm , the probability that it will end up at Rfi na i- 
Contours show 50% and 75% probability. The red dashed lines in the N300 panels show the radial extent (in scale lengths) of NGC 300 
studied in the current work. Yellow dashed lines show one-to-one correspondence between formation and final radii, i.e. the expected 
position with no migration. The N300 model shows less migration than the MW model. 



(|Sellwood fc Carlberd[l9 84), and the only way to cool a 
stellar disk is by star formation, which repopulates cir- 
cular orbits with young stars. Hence, if the majority of 
a disk is old, as in NGC 300, we can expect that it ex- 
perienced radial redistribution early in its evolution, but 
not in recent years. Because our time resolution at old 
times is poor, we would not be sensitive to such early 
evolution. 

In contrast to NGC 3 00, the majority of sta rs in M33 
formed more recently ( William s et alJ l2009a|) . so M33 
may have experienced more migration in recent times. 
M33's interaction with M31 may also have driven en- 
hanced spir al structure, which in turn drives more mi- 
gration (e.g. lQuillen et~aT1[2009l) . NGC 300 is in a much 
more isolated environment. These factors combined sug- 
gest that one would expect M33 to have experienced 
more migration than NGC 300, not less. 



In summary, while migration will tend to erase the sig- 
nature of disk growth, it is not the only way to obtain 
scale lengths which do not evolve as a function of age 
of stellar population. From our modeling, we can sur- 
mise that mass is an important factor in the extent of 
radial migration, with less massive systems experienc- 
ing less redistribution. Consequently, the lack of scale 
length evolution in NGC 300 is attributable to a factor 
other than radial migration, in this case most likely the 
rapid early growth of the disk. 

Further evidence that NGC 300 has not experi- 
enced much migration comes from the ki n emati cs of 
the globular cl uster systems. lOlsen et"al1 (|2004D and 
iNantais et al.l ()2008[ ) found that the globular clusters 
in NGC 300 had kinematics matching that of the Hi 
disk, indicating that their present location is where they 
formed. Since both heating and migration are caused by 
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Fig. 11. — Mean metallicities corresponding to SFHs in Figure l4l Errors are calculated as in Figure [4] 



spiral arms, this finding would suggest a lack of migration 
as well. In a more massive galaxy, NGC 253, asymmet- 
ric drift of the globular clusters indicates stronger radial 
diffusion. This result agrees with the prediction of the 
models described here in suggesting that more massive 
galaxies experience greater migration. The old globular 
clusters with disk kinematics in NGC 300 also fits with 
the interpretation that the majority of stars formed early. 

4.2.4. Comparison with M33's scale length evolution 

We considered whether the NGC 300 data is consistent 
with the stronger scale length evolution seen in M33, such 
that the observed weaker evolution is due solely to the 
shallower depth of the NGC 300 CMDs. In other words, 
if M33 were at the distanc e of NGC 300, would the results 
of lWilliams et al.l (|2009aD be recoverable? Unlike the test 
in cj4.2.1[ in which we saw that the SFH of M33 is not 
compatible with the NGC 300 data, here we are looking 
at the effect of depth on the M33 data. 

We calculated photometry of the M33 data using only 
one 1300s exposure in F &Q&W and FilAW for ea ch of the 
ACS fields analyzed in iWilliams et~aT1 (|2009aD . When 
deriving the SFH we considered only the portion of the 
CMD that was 2 magnitudes above the 50% completeness 
limit, to mimic the effect of placing M33 at the distance 
of N GC 300. We found a sm aller scale length evolution 
than [Williams et al.l (|2009af) . with the scale length in- 
creasing from 1.0 ± 0. 2 kpc at 8 Gyr to 1.3 ± 0.1 kpc at 
0.1 Gyr. In contrast, IWilliams et all ([2009a) found an 
increase of 1.0 ± 0.1 at 10 Gyr to 1.8 ± 0.1 at 0.6 Gyr 
with the deeper data. Thus, the reduced depth can lead 
to an underestimate of the scale length evolution; it may 
be possible that the evolution of NGC 300 is comparably 
dramatic to that seen in M33, but that we are not sen- 
sitive to it due to the shallower depth of the NGC 300 
data. On the other hand, the inferred evolution for the 
shallower M33 data is still somewhat larger than we ob- 



served in NGC 300. We emphasize that the difference in 
mean age of the two galaxies, as seen in H4.2.11 is still 
secure. 

4.3. Metallicity 

In most models of galaxy evolution, the buildup of a 
stellar population is accompanied by an increase in the 
mean stellar metallicity. Outside the Milky Way, metal- 
licity is often measured only for the youngest stellar pop- 
ulations, using either Hn regions or atmospheres of A 
and B stars (with the latter method being used only for 
the nearest galaxies). Here we discuss the present-day 
metallicity structure of NGC 300 and compare it to its 
past evolution as derived from stellar populations. 

4.3.1. Present Day Metallicity 

A radial metallicity gradient for NGC 300 in the 
gas and young s tars has been reported by a num- 
ber of authors (IPagel et al.l 1197! 
Edmunds fc Pagell 1198 
Zaritskv et al 



1983; 



Web ster fc Smith 
IDeharveng et al. 
Urbaneia et al.l 120051: 



1988 ; IZaritsky et all Tl994l: 
Kudritzki et al.l 120081: iBresolin et al.l 12009ft . All these 
studies found that metallicity is highest in the center 
and decreases toward larger radii, although the overall 
metallicity depends on the calibration method used. 
The results best suited for comparison with metallic- 
ities derived fr o m ou r ste llar populations are th ose of 
lUrbaneia et~aT1 (|2005l) and Kudritzki et al.l (|2008l ). since 
they determine metallicities using stellar spectroscopy of 
individual young A and B stars from the galactic center 
out to ~ 6.8 kpc, and are thus measuring the stellar, 
rather than gas-phase, metallicity grad i ent. T he best- fit 
abundance gradient of Kudritzki et all (|2008| ) is 



[Z] = (-0.06 ± 0.09) - (0.078 ± 0.021)r 



(2) 



where r is radius in kpc (assuming the distance used in 
this paper of 2.0 Mpc) and [Z] is an average metallicity 
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TABLE 2 



ANGST DATA: . 


AGES, SFRS, AND METALLICITIES FOR EACH RADIAL 








BIX 




Radius 


Area 


Age 


SFR 


[M/H] 


(kpc) 


(kpc 2 ) 


(Gyr) 


(M yr" 1 kpc" 2 ) 




0-0.9 


1.40 


0.004-0.079 


0.0110 ±0.0013 


-0.24 ±0.05 






0.079-0.13 


0.0069 ± 0.0009 


-0.24 ±0.05 






0.13-0.25 


0.0054 ± 0.0011 


—0.25 ± 0.05 






0.25-0.40 


0.0511 ± 0.0109 


-0.23 ± 0.05 






0.40-10 


0.0219 ± 0.0085 


-0.51 ± 0.05 






10-14 


0.1098 ± 0.0291 


-0.92 ± 0.10 


0.9-1.8 


2.53 


0.004-0.079 


0.0189 ±0.0023 


-0.34 ±0.05 






0.079-0.13 


0.0117 ±0.0018 


-0.34 ±0.05 






0.13-0.20 


0.0100 ± 0.0009 


-0.34 ±0.05 






0.20-0.32 


0.0102 ±0.0012 


-0.35 ±0.05 






0.32-0.50 


0.0410 ± 0.0107 


—0.34 ± 0.05 






0.50-10 


0.0041 ± 0.0036 


-0.47 ±0.06 






10-14 


0.1358 ± 0.0488 


-0.78 ± 0.16 


1.8-2.7 


1.98 


0.004-0.079 


0.0089 ± 0.0015 


-0.38 ± 0.05 






0.079-0.13 


0.0069 ± 0.0016 


-0.38 ±0.05 






0.13-0.20 


0.0077 ± 0.0008 


-0.38 ±0.05 






0.20-0.32 


0.0081 ± 0.0005 


-0.38 ± 0.05 






0.32-0.50 


0.0162 ± 0.0019 


-0.38 ± 0.06 






0.50-1.6 


0.0064 ± 0.0009 


—0.38 ± 0.06 






1.6-10 


0.0046 ± 0.0043 


-0.51 ± 0.08 






10-14 


0.0452 ± 0.0141 


-0.82 ± 0.14 


2.7-3.6 


1.90 


0.004-0.079 


0.0124 ± 0.0024 


-0.47 ±0.11 






0.079-0.13 


0.0074 ± 0.0008 


-0.48 ±0.11 






0.13-0.20 


0.0054 ± 0.0010 


-0.48 ±0.11 






0.20-0.40 


0.0045 ± 0.0005 


-0.48 ±0.11 






0.40-0.63 


0.0090 ± 0.0012 


-0.48 ±0.11 






0.63-2.0 


0.0049 ± 0.0025 


—0.49 ± 0.12 






2.0-10 


0.0035 ± 0.0029 


-0.60 ± 0.06 






10-14 


0.0162 ± 0.0060 


-0.92 ± 0.16 


3.6-4.5 


1.83 


0.004-0.100 


0.0050 ± 0.0009 


-0.41 ±0.07 






0.100-0.158 


0.0034 ± 0.0004 


-0.42 ±0.07 






0.158-0.25 


0.0035 ± 0.0004 


-0.42 ±0.07 






0.25-0.40 


0.0044 ± 0.0008 


-0.41 ±0.07 






0.40-0.79 


0.0045 ± 0.0005 


-0.43 ± 0.07 






0.79-3.2 


0.0022 ± 0.0011 


—0.43 ± 0.07 






3.2-10 


0.0042 ± 0.0031 


-0.68 ±0.04 






10-14 


0.0052 ± 0.0025 


-1.29 ±0.31 


4.5-5.4 


1.35 


0.004-0.079 


0.0027 ± 0.0006 


-0.58 ±0.10 






0.079-0.158 


0.0015 ±0.0002 


-0.58 ±0.09 






0.158-0.25 


0.0017 ± 0.0003 


-0.58 ±0.10 






0.25-0.40 


0.0014 ± 0.0002 


-0.58 ±0.09 






0.40-0.79 


0.0022 ± 0.0005 


-0.59 ±0.10 






0.79-1.6 


0.0033 ± 0.0003 


-0.60 ±0.09 






1.6-5.0 


0.0023 ± 0.0012 


-0.71 ±0.07 






5.0-14 


0.0011 ±0.0005 


-1.21 ±0.23 



TABLE 3 

Fits to the surface density and metallicity gradient 



Age 


So 


Rd 


Zo 


AZ/Ar 


(Gyr) 


(M pc- 2 ) 


(kpc) 




(kpc- 1 ) 



0.004 506.946 ± 106.331 1.21 ±0.10 -0.23 ± 0.05 -0.072 ± 0.019 

0.1 507.602 ± 106.786 1.21 ±0.10 -0.23 ± 0.05 -0.075 ± 0.019 

1.0 513.618 ± 110.948 1.16 ±0.10 -0.40 ± 0.05 -0.057 ± 0.018 

5.0 586.142 ± 129.738 0.95 ± 0.09 -0.52 ± 0.05 -0.078 ± 0.020 

10 443.474 ± 110.079 0.97 ±0.10 -0.89 ±0.11 -0.059 ± 0.049 



based on a combination of multiple abundance measure- 
ments (mostly Fe, Ti, and Cr lines). We adopt Equa- 
tion [5] as the best estimate of the present-day metallicity 
gradient. The Hii-region metallicit ies from temperature - 
sensitive emission lines presented in lBresolin et al.l (|2009[ ) 
also agree remarkably well with these stellar metallicities. 

4.3.2. Metallicity Evolution from HST Data 



The metallicit y grad ient measured in young stars by 
iKudritzki et~aT1 J2008) can be compared to the present- 
day metallicities inferred from the observed stellar pop- 
ulations. The CMDs contain information on the metal- 
licities as well as the ages of stars, espec ially for stars 
off the main sequence (jGallart et al.l 120051 ). A full suite 
of isochrones with varying age and metallicity is used to 
fit each CMD, so that in the final SFH each time bin is 
associated with a particular metallicity. The full metal- 
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r (kpc) 

Fig. 12. — Solid black li ne: best-fit metallicity gradient obtained 
from young massive stars (Kudritzki ct al. 2008), with shaded re- 
gion showing uncertainties. Blue circles are metallicities measured 
from our CMDs in each radial bin at the present day (i.e., the 
mean metallicity for the past 100 Myr); green, orange, and red 
circles are metallicities for the 1—5 Gyr, 5— 10 Gyr, and 10-14 Gyr 
bins, respectively. Solid lines of corresponding colors show the best 
linear fit for each set of points. Errors are calculated as in Figure|4] 
Mean metallicities for the entire stellar population are shown with 
black triangles, with a dashed black line showing the best linear fit. 
Points have been offset in radius to avoid overlapping error bars. 



licity history for all radial bins is shown in Figure 111! 
Uncertainties in the metallicity may result from several 
factors: 1) an age-metallicity degeneracy in the RGB, 2) 
photometric errors in the CMD, 3) uncertainties in the 
reddening, 4) uncertainties in modeling stars on the in- 
stability strip, and 5) errors in the isochrones, especially 
for supergiants and AGB stars. We address these uncer- 
tainties by constraining the metallicity not to decrease 
significantly with time, which eliminates unphysical so- 
lutions to the SFH. 

We tested the effect of this parameter by deriving the 
SFH without it as well. Not constraining the metallicity 
results in lower [M/H] values for intermediate age stars 
(1-10 Gyr), especially in the most crowded regions, but 
the effect on the SFR is minimal. The trends in SFH 
reported in this paper do not depend on whether or not 
the metallicity constraint is enforced. 

The derived SFH from the ANGST data reproduces a 
metallicity gradient in the present-day m etallicity (4-100 
Myr) time bin. In Figure[T2]we show the lKudritzki et all 
(2008) metallicity gradient along with our derived metal- 
licity gradient. Fit parameters for the linear fits are given 
in Table |3] The MATCH-derived present-day metallic- 
ities are generally consistent with the observed values 
within the error bars, although the metallicities are all 
systematically low and fall near the minimum value in 
the observed uncertainty range. The mean metallicities 
of the entire population are lower as well, as would be 
expected for a stellar population that built up gradually 
over ~ 10 Gyr. However, the mean stellar metallicity 
still shows a radial gradient at all epochs. 

Since the color of the main sequence is not significantly 
affected by metallicity, our CMD-fitting code may not be 
sensitive to a metallicity increase at very recent times. 
The metallicities we find are therefore a lower limit on 
the current values (at least in the central regions where 
the infall of unenriched gas is assumed to have ended). 



This insensitivity is especially true in the innermost ra- 
dial bins, in which the low current SFRs and the effects 
of crowding make the present-day metallicity especially 
difficult to determine. Not coin cidentally, these values 
fall more sig nificantly below the IKudritzki et all (2008) 
results than do the values for the outer disk. 

The change in the metallicity gradient over time can 
be seen in the mean metallicity derived during multiple 
time bins, also shown in Figure 1121 Stellar population 
ages of 4-100 Myr, 1-5 Gyr, 5-10 Gyr, and 10-14 Gyr 
are shown. (The metallicity values for age 100 Myr - 1 
Gyr are essentially identical to the present-day values.) 
Results of linear fits to the gradients are presented in 
Table [3] Outside 2 kpc, the slope of the metallicity gra- 
dient has remained fairly constant over time, with the 
overall metallicity increasing. In the central 2 kpc, we 
see a significant metallicity increase from the 10-14 Gyr 
bin to the 5-10 Gyr bin, which corresponds to the epoch 
when significant star formation activity was taking place 
in the galactic center. 

4.3.3. Comparison with M33's metallicity 

Figure shows that NGC 300 has both a signifi- 
cant radial metallicity gradient and a significant increase 
in metallicity with time. Although we constrained the 
metallicity not to decrease with time, there was no bar- 
rier to its remaining constant in deriving the SFHs, or 
to its taking arbitrary slopes with radius. In contrast, 
M33 has a shallower gradient and metallicity which has 
remained nearly constant with time (jMagrini et al.|[2009l 
Holtzman et al. 2009, in prep.). 

These results are consistent with chemical evolution 
models. Ongoing star formation in M33 may indicate 
gas infall, which can dilute the metallicity, keeping it con- 
stant with time despite the continuing enrichment from 
star formation. In contrast, NGC 300 formed a small 
fraction of its stars at late times, which may indicate a 
less important role for ongoing gas accretion, and thus 
chemical evolution that behaves more like a closed box. 
This may also be a counterexample to the "downsizing" 
model, since NGC 300 is a fairly low-mass disk galaxy. 

5. CONCLUSIONS 

We have presented resolved stellar photometry of 3 
HST/ ACS fields in a continuous radial strip from the 
center of NGC 300 to 5.4 kpc, as well as additional 
archival fields scattered throughout the disk. For both 
the ANGST (continuous) and archival fields, we have di- 
vided the stars into 6 radial bins of width 0.9 kpc and 
used CMD-fitting to derive the SFH in each bin. We can 
resolve the red clump in the majority of the radial bins. 
We find that the disk of NGC 300 has a dominant old 
population throughout the observed region: looking at 
the disk as a whole, ~80% of stars are older than 6 Gyr. 
However, the outer parts of the disk have a higher per- 
centage of young stars than the inner parts, consistent 
with inside-out growth. In the inner regions, > 90% of 
the stars arc older than 6 Gyr, while in the outermost ra- 
dial bin (4.5 < 5.4 kpc), only ~40% of stars are this old. 
Comparison of formation and final positions of stars in an 
TV-body simulation of an NGC 300-like galaxy indicates 
that the effects of migration should not be substantial in 
a galaxy of this size. 
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We calculated the surface density profile of the disk at 
each time step by summing over the stellar mass formed 
and fitting the profile to derive the corresponding scale 
lengths. According to inside-out growth, the scale length 
of the disk should increase with time. While we found 
this to be the case, the increase is relatively small, from 
~ 1.1 to ~ 1.3 kpc over the lifetime of the disk. In 
contrast, the scale lengths of M33 found with an identical 
method show a much greater i ncrease, from 1.0 kpc to 
1.8 kpc (|Williams et alj I2009al) . more in line with the 
theoretical predictions of I Mo et al.l (|1998f) . 

Although M33 is nearly a twin to NGC 300 in Hub- 
ble type and mass, there are differences between the two 
galaxies which may be related to their different SFHs. 
The primary difference in morphology is the presence of 
a disk break in M33. Breaks have been shown to be 
extre mely common in disk galaxies (jPohlen fc Truiiliol 
2006), and M33 is no exception, with a break at 8 kpc 



dFerguson et al.l 20071) . The scale length increase seen by 



2009a|) only holds inside the disk break; 



iWilliams et al.l 
outside the break the opposite trend is seen, with a de- 
creasing scale length with time. In contrast, NGC 300 is 
unusual in that it has no disk br eak out to at least 14 kpc 
Bland-Hawthorn et al. 2005). Sanchcz-Blazqucz et al.l 



2009) suggest that pure exponential profiles may be a 



feature of galaxies that have undergone less radial mix- 
ing, which is consistent with our analysis in £14.2.31 that 
NGC 300 has undergone relatively little recent migration. 

Environmental factors may also be important, 
as NGC 300 is isolated from other large galax- 
ies while M33 appears to be interacting with M31 
(iBraun fc Thilkeril2004; lBekkill2008l iMcConnachie et al.l 
120091: iPutman et all 12009(1 . An influx of gas onto M33 
in recent times may have triggered star formation in the 
outskirts of M33 and contributed to the growth of the 
disk. 

Finally, we note that despite the unbroken exponen- 
tial profile of NGC 300, it is not without features in the 
very outer d i sk (b eyond the extent of our observations) : 
IVlaiic et al.l ((2009) find a change in the abundance gradi- 
ent at ~ 10 kpc. They propose two possible explanations 
for the upturn in the gradient: radial mixing, which we 
conclude is unlikely to be the sole effect based on the 
simulation presented in H4.2.31 and an accretion scenario 



in which the outer regions of the disk form stars later 
and become enriched, flattening the gradient. We sug- 
gest that additional HST observations of resolved stars 
farther out in NGC 300 would provide useful constraints 
on the total disk scale length and evolution. 

From our derived SFH for NGC 300, we find a present- 
day metallicity gr adient roughly consistent with ob- 
serva tional results (jKudritzki et alj|2008t IBresolin et al.l 
2009). Additionally, we find that the metallicity has in- 
creased with time in all radial bins, suggesting a lack of 
infall of unenriched gas. This is broadly consistent with 
our finding that the majority of stars in NGC 300 formed 
prior to 6 Gyr ago, as more recent gas infall would likely 
have triggered more enhanced recent star formation. As 
M33 has been shown to have little evolution in metal- 
licity and more of its stars formed recently, it may have 
experienced more gas infall at later times. Despite the 
visual similarities between these two galaxies, they seem 
to have markedly different histories. 
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APPENDIX 



ARCHIVAL DATA 

In this section we present the results of analyzing archival fields in NGC 300 (see H2.1j) in the same manner as the 
ANGST data. CMDs for the archival data are shown in Figure [131 The quality cuts used for selecting stars were 
identical to those described for the ANGST data in £12.21 Although the archival fields are not radially aligned, the 
stars in these fields can be sorted into the same 6 bins as the new observations. These fields are not evenly spaced in 
radius, as shown in Figure [TJ and were placed on regions of high current star formation. We include only stars up to 
a galactocentric distance of 5.4 kpc, the radial extent of the ANGST data. 

As with the ANGST data, we fixed the distance for these SFH derivations at m — M = 26.43. The results of the 
fits were combined into time bins identical to those used for the ANGST data. For the archival data, mean extinction 
values for each bin are, from inner to outer: A v = 0.30±0.05, 0.20±0.05, 0.10±0.05, 0.10±0.05, 0.10±0.05, 0.13±0.07. 
The higher extinction values for the archival data are likely due to the location of these fields in regions of higher star 
forma tion, which should correspond with increased dust content that affects older stars as well. Indeed, IRoussel et al.l 
(2005) studied extinction in NGC 300 and found that extinction is variable for young clusters. 

Reassuringly, the SFH from the archival data, shown in Figures Q3] and [131 is very similar to that derived from the 
ANGST data, despite the different filters and depths. Cumulative SFHs for archival data divided into radial bins are 
shown in Figure [TBI and integrated for the whole galaxy in Figure [TTj Metallicities are shown in Figure [18] Table [4] 
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Fig. 13. — CMDs of radial bins from archival data. Quality cuts are identical to those described in FigureO 
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Fig. 14. — Lifetime SFHs of radial bins from archival data (solid lines). Errors are calculated as in Figure[4] For comparison, the ANGST 
SFHs arc shown as dashed lines, with uncertainties shown as shaded regions. 
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Fig. 15. — Recent (< 1 Gyr) SFHs of radial bins from archival data (solid lines). Errors are calculated as in Figure [J] For comparison, 
the ANGST SFHs are shown as dashed lines, with uncertainties shown as shaded regions. 
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Fig. 16. — Cumulative star formation for all radial bins for archival data. Bin edges have been offset to avoid overlapping error bars. 



gives the full SFH for the archival data. 
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TABLE 4 

Archival data: ages, SFRs, and metallicities for. each radial 

bin 
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